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Abstract One of the main ingredients in most quantum information protocols is a reliable 
source of two entangled systems. Such systems have been generated experimen- 
tally several years ago for light [1-5] but has only in the past few years been 
demonstrated for atomic systems [6-9]. None of these approaches however in- 
volve two atomic systems situated in separate environments. This is necessary 
for the creation of entanglement over arbitrary distances which is required for 
many quantum information protocols such as atomic teleportation [10, 11]. We 
present an experimental realization of such distant entanglement based on an 
adaptation of the entanglement of macroscopic gas samples containing about 
lO^'^ cesium atoms shown in [8]. The entanglement is generated via the off- 
resonant Kerr interaction between the atomic samples and a pulse of light. The 
achieved entanglement distance is 0.35m but can be scaled arbitrarily. The fea- 
sibility of an implementation of various quantum information protocols using 
macroscopic samples of atoms has therefore been greatly increased. We also 
present a theoretical modeling in terms of canonical position and momentum 
operators X and P describing the entanglement generation and verification in 
presence of decoherence mechanisms. 



1. Introduction 

Ever since Einstein, Podolslcy, and Rosen in their seminal paper from 1935 
[12] introduced the possibility of entangling two quantum system, entangle- 
ment has been viewed as one the most curious and spectacular phenomena in 
quantum mechanics. In the past few years the role of entanglement in quantum 
mechanics has shifted dramatically from being a fundamental test of the foun- 
dation of the entire quantum mechanical theory to being a technical resource 
in the rapidly developing field of quantum information. Thus, the hunt is on 
for reliable sources of entanglement. These have been available for discrete 
as well as continuous states of the electromagnetic field. However, entangled 
states of material particles have presented a greater experimental challenge. 
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Macroscopic samples of atoms as a resource of entanglement have attracted 
a lot of attention in recent years because of their relatively simple experimental 
realization (works at room temperature) and robustness to single particle de- 
coherence. In [8] entanglement of this kind was accomplished. However, the 
two samples were located only 1cm apart in the same shielded environment. 
This meant that this implementation did not incorporate the important feature 
in e.g. the teleportation protocol, that the distance of teleportation given by the 
separation of the entangled systems could be arbitrary. We have created two 
separate environments each containing a gas sample of cesium atoms at room 
temperature. As we will see, we have successfully created entangled states 
between two such systems being 0.35m apart. In the present paper we will 
also focus on a better method to verify the generation of entangled states as 
compared to the experiment in [8]. 

2. Light Atom Interaction 

In this section we introduce the physical systems involved in the experiment, 
i.e. we introduce the atomic spin samples and the polarization state of laser 
pulses interacting with each other Based on the equations of motion we will in 
Sec. 3 explain how the interaction can be utilized for entanglement generation. 

Atomic System 

Our atomic system is composed of two separate samples of spin polarized 

cesium vapour placed in paraffin coated glass cells at room temperature. Ce- 
sium has a hyperfine spht ground state with total angular momentum F = 3 
and F = 4, the latter being our atomic quantum system of interest. Having 
a macroscopic ensemble of atoms (around 10^^) we will define the theoreti- 
cally discrete but effectively continuous collective spin variables Jk = J2ijk' 
where k = x,y,z and i denotes the individual atom. These will retain regu- 
lar angular momentum commutation relations, [Jy, Jz] = iJx- The Heisenberg 

uncertainty relation will then lead to Var(Jy)Var(J2) > i| ^-^x^ P- We will be 
interested in the state in which practically all atoms are in the F = 4:, mp = A 
state with x as quantization axis. All the relevant interactions only involve 
minute changes in the macroscopic spin orientation. can then be regarded 
as a constant classical number. If all atoms are independent this will give rise 
to a minimum uncertainty state called the coherent spin state (CSS) with: 

Var(4) = Var(J,) = = 2iVatoms. (CSS) (1) 

We therefore see that the variance, often referred to as the projection noise, of 
the CSS will grow proportionally to the number of atoms. This scaling is man- 
ifestly quantum and will be important for the verification of the entanglement. 
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Light System 

In complete analogy to the atomic collective spin variables, the polarization 
state of a pulse of light can be described by a vector, the so-called Stokes 
vector. For light propagating along the z-axis we define 

Sx{t) = ^{^x- (t>y) , Sy{t) = ^ (045 - 0135) , Sz{t) = ^ (0+ - 0_) , 

(2) 

where 4>x, 4>y the photon fluxes of x and y-polarized photons, 045, ^135 are 
photon fluxes measured in a basis rotated 45° with respect to the x, y-axes, and 
0+, 0_ refer to photons with a+ and (T_ -polarization. In our experiments the 
light will with very good approximation be linearly polarized along the x-axis. 
Then Sx{t) can be described by a classical c-number. The Sy{t) and Sz{t) 
operators will contain the interesting quantum variables. The Stokes opera- 
tors defined above have dimension time~^. This is convenient for describing 
light/matter interactions as we will see below. But for entanglement generation 
or quantum information protocols in general it is more convenient to consider 
entire pulses which are time integrated versions of the above. It can be shown 
that 

/ Sy{t)dt, [ Sz{t)dt\ = i [ Sx{t)dt = nph/2, (3) 
Jo Jo J Jo 

where the last equality holds in our case for strong hnear polarization along 
the X-axis. From this commutator we derive the variance of the minimum 
uncertainty state called the coherent state: 

Var Sy{t)d?j = Yai Sz{t)d?j = Up^/A. (shot noise) (4) 

Note again the characteristic linear quantum scaling with the number of parti- 
cles nph. We often refer to the variance of the coherent Ught state as the shot 
noise level. 

For completeness we will note that for coherent light states it can be shown 
[13] that 

Sy{t)Sy{t')) = {Sz{t)Sz{t')) = ^6{t - t'). (5) 

The usefulness of this will be shown below. In our experiment we will use 
Sy detections to create entanglement between atomic samples. If initially x- 
polarized Ught is rotated a small angle 9 around the axis of propagation (z) we 
get (^Sy'j = 2Sx9. This is why an Sy detection is sometimes referred to as a 
polarization rotation measurement. 



4 



Interaction 

We couple our light and atomic system by tuning a laser beam off -resonantly 
to the 6S'i/2 6P3/2 dipole transition in cesium. This leads to the following 
equations of interaction: 

ST(t) = Sl;^{t)+aSMt), (6) 
= STit), (7) 

^Jy{t) = ajJTit), (8) 
§-^Ut) = 0, (9) 



where a = — g^xA- ^ beam cross section, A is the detuning (red posi- 
tive), A is the optical wavelength, and 7 is the natural linewidth of the excited 
state QPs/2- In and out refer to light before and after passing the atomic sam- 
ple, respectively. The above equations have been derived carefully in [13] but 
we will give a short physical explanation here. 

First of all, the interaction is refractive in nature (the absorption of the off- 
resonant light is negligible). It is convenient to consider the incoming linearly 
polarized light in the (T+ and a- basis. The phase shift of a and a a- photon 
propagating through atoms will be different if there is a spin component 
along the propagation direction. For instance (quantized along z) an atom in 
the mp = F magnetic sub-level couples strongly to cr+ photons and weakly 
to cr_ photons. For the mp = —F sub-level the situation is reversed. The 
differential phase shift of a+, a- photons turns out to depend linearly on 
which leads to a polarization rotation of the incoming linearly polarized Ught 
proportionally to (also known as Faraday rotation). Eq. (6) is a first order 
approximation of this effect. 

The different coupling strengths for different sub-levels also lead to a Stark 
shift of atomic levels depending on the sub-level quantum number mp and 
the incoming light polarization. Integrated over time, the Stark shifts lead 
to different phase changes of the magnetic sub-levels which changes the spin 
state. If for instance there are more (7_|_ than a- photons, the mj? = F sub- 
state will be affected more than the = —F state. The amount of a-\- and 
cr_ photons is measured by the 5"^ operator. It turns out that the spin state 
evolution can be described as a rotation of the spin around the z-axis by an 
amount proportional to Sz- This is given to first order in Eq. (8) (the rotation 
is so small that Jx is unaffected). 

In the interaction process the cr_|_ and (T_ photons experience phase shifts 
but are not absorbed. In the off-resonant limit the flux of a-^- and <t_ photons 
are individually conserved leading to Eq. (7). By conservation of angular mo- 
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mentum along the z-direction this leads to the constancy of Jz expressed by 
Eq. (9). 

We see from Eqs. (6) and (9) that in the case of a large interaction strength 
(i.e. if aSxJz dominates S"™) a measurement on 5",°"* amounts to a measure- 
ment of Jz without destroying the state of Jz- This is termed a Quantum Non 
Demolition (QND) measurement of Jz- Using off-resonant light for QND mea- 
surements of spins has also been discussed in [14, 15]. We note that Eq. (8) 
implies that a part of the state of light is also mapped onto the atoms. This 
opens up the possibility of using this sort of system for quantum memory. One 
step in this direction is discussed in [13, 16]. 

Adding a Magnetic Field 

In the experiment a constant and homogeneous magnetic field is added in the 
x-direction. We discuss the experimental reason for this below. For our model- 
ing the magnetic field adds a term Hb = ^Jx to the Hamiltonian. This makes 
the transverse spin components precess at the Larmor frequency Q depending 
on the strength of the field. Introducing the rotating frame coordinates: 

/ j; \ _ / cosm sm{nt) \(Jy\ 

WW ~ V -^Hm cos(j^t) j\jz J ^ ^ 

we can easily show that Eqs. (6)-(9) will transform into: 

S'^\t) + a5a; (^J^{t) sm{nt) + J'z{t) cos(Qi)) , (11) 
ST{t), (12) 
aJ^S'^{t) cos{m), (13) 

aJxS'^{t)siii{9.t). (14) 

Thus, the atomic imprint on the light is encoded in the Jl-sideband instead of 
at the carrier frequency. The advantage of this added feature is threefold. The 
first and perhaps most important advantage is that lasers are generally a lot 
more quiet at high sideband frequencies compared to the carrier. A measure- 
ment without a magnetic field will be a DC measurement and the technical 
noise would dominate the subtle quantum signal. Secondly, the 5-field in- 
troduces a Larmor splitting of the magnetic sublevels of the hyperfine ground 
state multiplet, thus lifting the degeneracy. This will introduce an energy bar- 
rier strongly suppressing spin flipping collision. The lifetime of the atomic 
spin state is consequently greatly increased. The last advantage is that as long 
as the measurement time is longer than l/fi Eq. (1 1) enables us to access both 
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Jy and at the same time. We are of course not allowed to perform non- 
destructive measurements on these two operators simultaneously since they 
are non-commuting. This is also reflected by the fact that neither Jy nor Jz 
are constant in Eqs. (13) and (14). In section 3 we shall consider two atomic 
samples and the third advantage becomes evident. 

3. Entanglement Creation and Modeling 

In this section we define what is meant by entangled states of atomic samples 
and we adapt the equations of motion from previous sections for this purpose. 
We will derive a simple model for the entanglement creation and describe how 
to verify that the states created really are entangled. In [17] a much more 
elaborate description of entanglement creation is given. 

Entanglement Criterion 

Let us here state the criterion to fulfil in order to prove the generation of 

entangled states. Since entanglement is the nonlocal interconnection of two 
systems we need to have two atomic samples A and B. Entanglement is usually 
defined in terms of density matrices so that A and B are entangled if they are 
connected in such a way that it is impossible to write the total density matrix 
as a product, ptot = J^PiPAiPBi- For our continuous variable system we have 
the experimentally practical criterion derived from the above definition in [18]: 



Var( J^i + J^s) + Var(J^i + J',^) < 2 (15) 



where we have assumed both samples to be macroscopically oriented with 
same magnitude Jx- The two samples are indexed by 1 and 2. Comparing to 
Eq. (1) we get an equality for two independent atomic samples in the CSS. 
To have entanglement we thus need to know the sums of the spin compo- 
nents along the y- and the z-directions better than we could ever know each 
of the spin projections by itself. It is now interesting to examine the commu- 
tator between the sums, -\- Jy2i J'zi ~^ z2 — '^{.Jxi "I" Jx2)- A non-zero 
commutator means that increasing our knowledge of one component will au- 
tomatically decrease our knowledge of the other, thus making our attempts to 
break the inequahty of Eq. (15) futile. Now comes the trick that makes en- 
tanglement generation possible in our experiment. Assume Jxi and Jx2 to be 
equal in magnitude but opposite in direction. The commutator will then be- 
come zero and we can at least theoretically measure both components with 
arbitrary precision, thereby satisfying the entanglement criterion of Eq. (15). 
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Two Oppositely Oriented Spins 

Inspired by the above we will from now on assume Jxi = —Jx2 = Jx- 
We will re-express the equations of motion (1 1)-(14) for two samples in a way 
which is much more convenient for the understanding of our entanglement 
creation and verification procedure. We introduce position and momentum Uke 
operators X and Pto describe pulses of light and the atomic systems. This is a 
more abstract but hopefully also more well known and intuitive way to express 
the interactions creating the entangled states. 

For two atomic samples we write equations of motion: 

S'^it) + aSx + j;2(i)] sin(Oi) (16) 

a{Jxi + Jx2)S'^\t) cos{nt) = 0, (17) 
a{Jxi + Jx2)S'^{t) sin(m) = 0. (18) 



^{J'ylit) + J',^{t)) = 

^AJ'At) + J'At)) = 



The fact that the sums J'yiit) + Jy2if) and J'^iit) + <^z2(*) ^^^^ ^^''^ tinie 
derivative relies on the assumption of opposite spins of equal magnitude. The 
constancy of these terms together with Eq. (16) allows us to perform QND 
measurements on the two sums. We note that each of the sums can be accessed 
by considering the two operators 



cos{^t)dt 



^ S'^cos{i}t)dt + + J'At)), (19) 



Jo 

/ 57*sin(m)dt= / S'^sm{nt)dt + ^{JL{t) + JUt))- (20) 
Jo Jo 2 



We have used the fact that cos'^{ilt)dt sm'^{i}t)dt « 1/2 and that 

Jq cos{flt) sm{Q,t)dt « 0. Each of the operators on the left hand side can be 
measured simultaneously by making a S'j^-measurement and multiplying the 
photocurrent by cos{flt) or sin(nt) followed by integration over the duration 
T. The possibihty to gain information about J^i(t) + Jy2{t) "^'zii^) + 
Jz2{t) enables us to break the inequality (15). At the same time we must 
loose information about some other physical variable. This is indeed true, the 
conjugate variables to these sums are J^2(*) ~ j'zii^) ~ •^«2(0' 



respectively. These have the time evolution 

d 



^(i^iW - J'At)) = 2ajJT{t)sm{nt). 



(21) 
(22) 



We see how noise from the input -variable is piling up in the difference 
components while we are allowed to learn about the sum components via Sy 
measurements. The above equations clearly describe the physical ingredients 
in play but the notation is cumbersome. Therefore we define new operators. 
For the atomic system we take 



Xai 
Pai 

Xa2 
Pa2 



T' ~ V 
V + V 
V - V 

\/2 Jj; 



New Ught operators will be 

Pli = 

Xl2 = 
Pl2 = 



Jo 



(t) cos{nt)dt, 



Sz{t) cos{^lt)dt, 



Sy{t) sm{nt)dt, 



Sz{t) sm{m)dt. 



(23a) 
(23b) 
(23c) 
(23d) 

(24a) 
(24b) 
(24c) 
(24d) 



Each pair of X, P operators satisfy the usual commutation relation, e.g. we 



have 



Xli,Pli 



i. All previous equations now translate into 

y-out _ vin , ^ pin 



P'lOUt 
Li 



Lit 



-^Ai — ^Ai + I^^Li^ 

riout Din 

^Ai — ^Av 



(25a) 
(25b) 
(25c) 
(25d) 
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where we remember i = 1, 2 refer to the definitions above and not the two 
samples. The parameter describing the strength of light/matter-interactions is 
given by K = a\/ JxSxT. The limit to strong coupling is around k ^ 1. Note, 
we have two decoupled sets of interacting light and atomic operators. 

In the transition from Stokes operators to canonical variables in Eqs. (24a-d) 
the result (5) is a convenient tool for calculating variances. If for instance the 
input light state is the coherent vacuum state we have 



Sy{t)S^'{t'))cos{nt) cos{Q,t')dtdt' 



1 



'JxJ- JO JO ^ ' ^ 

(26) 

which is as expected. Likewise, if the two atomic samples are each in the 
coherent state we will derive e.g. Nw{Xa\) = 1/2. Coherent states of atomic 
or light systems as defined above correspond to what is known as coherent 
states of the X, P-operators. 

The entanglement criterion (15) written in X, P-language is 

Var(PAi) + Var(PA2) < 1. (27) 

We see that entanglement of the two atomic samples can be considered as so- 
called two mode squeezing. The uncertainty in the uncoupled pair of operators 
Pai and Pa2 is reduced on the expense of the increased noise in the operators 

Xai and Xa2- 



Entanglement Generation and Verification 

Now we turn to the actual understanding of entanglement generation and 
verification. Experimentally we perform the following steps (more details will 
be given in Sec. 4). First the atoms are prepared in the oppositely oriented 
coherent states corresponding to creating the vacuum states of the two modes 
{Xai, Pai) and {Xa2: Pa2)- Next a pulse of light called the entangling pulse 
is sent through atoms and we measure the two operators XJ^ and Xf^ with 
outcomes Ai and B\, respectively. These results bear information about the 
atomic operators Pa\ and Pa2 and hence we reduce variances Var(P^i) and 
Var(PA2)- To prove we have an entangled state we must confirm that the 
variances of Pai and Pa2 fulfil the criterion (27). That is we need to know the 
mean values of Pai and Pa2 with a total precision better than unity. For this 
demonstration we send a second verifying pulse through the atomic samples 
again measuring X£"* and X'^ with outcomes A2 and B2. Now it is a matter 
of comparing Ai with A2 and Bi with B2. If the results are sufficiently close 
the state created by the first pulse was entangled. 
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Now let us be more quantitative. The interaction (25a) mapping the atomic 
operators Pai out on Ught is very useful for a strong k and useless if «; ^ 1. We 
will describe in detail the role of k for all values. To this end we first describe 
the natural way to determine k, experimentally. If we repeatedly perform the 
first two steps of the measurement cycle, i.e. prepare coherent states of the 
atomic spins and performing the first measurement pulse with outcomes Ai 
and Bi, we may deduce the statistical properties of the measurement outcomes. 
Theoretically we expect from (25 a) 

{Ai) = {Bi)=0 and Var(^i) = Var(Bi) = ^ + y . (28) 

The first term in the variances is the shot noise (SN) of light. This can be 
measured in absence of the interaction where k = 0. The quantum nature of 
the shot noise level is confirmed by checking the linear scaling with photon 
number of the pulse, see Eq. (4). The second term arises from the projection 
noise (PN) of atoms. Hence, we may calibrate to be the ratio = PN/SN 
of atomic projection noise to shot noise of light. Theoretically has the 
linear scaling = aJxSxT with the macroscopic spin size Jx which must be 
confirmed in the experiment. 

Next we describe how to deduce the statistical properties of the state created 
by the entangling pulse. Based on the measurement results Ai and Bi of this 
pulse we must predict the mean value of the second measurement outcome. If 
K oo we ought to trust the first measurement completely since the initial 
noise of X™j is negligible, i.e. (^2) = Ai and (^2) = Bi. On the other 
hand, if k = we know that atoms must still be in the vacuum state such 
that (^2) = (-B2) = 0. It is natural to take in general (^42) = aAi and 
{B2) = aBi. We need not know a theoretical value for a. The actual value 
can be deduced from the data. If we repeat the measurement cycle N times 
with outcomes A^^\ b[^\ ^2 \ and .62*^ the correct a is found by minimizing 
the conditional variance 

Var(A2|^i)+Var(52|5i) = 



N 

i 

mm 



In order to deduce whether we fulfil the entanglement criterion (27) we com- 
pare the above to our expectation from (25a). For the verifying pulse we get 



P"ient / Tjent 
Ai ~ \ ^Ai 



= l + K2Var(Plf), 



(30) 
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where X]^^^'^'^'^ refers to the incoming light of the verifying pulse which has 
zero mean. P^^* refers to the atoms after being entangled. We see that the 
practical entanglement criterion becomes 



\ai{A2\Ai) + YaT{B2\Bi) = 1 + (Var(P|f ) + Var(Pl';^')) 
< 1 + = Var(Ai) + Var(Bi). 



(31) 



In plain English, we must predict the outcomes A2 and B2 with a precision bet- 
ter than the statistical spreading of the outcomes Ai and Bi with the additional 
constraint that Ai and Bi are outcomes of quantum noise limited measure- 
ments. 

Theoretical Entanglement Modeling 

Above we described the experimental procedure for generating and verify- 
ing the entangled states. Here we present a simple way to derive what we ex- 
pect for the mean values (i.e. the a-parameter) and for the variances Yai{P^l^). 

We calculate directly the expected conditional variance of A2 based on Ai : 



/ ( -v'Out,2nd -0-out,lst\ ^ \ 



(32) 



In the second step we assumed that the measurement is perfectly QND and 
without any decoherence, i.e. = P™^ . By taking the derivative with re- 
spect to a we obtain the theoretical minimum 



k2 



1 + ^2 

1 



(33) 



Var(A2|Ai) + \ax{B2\Bi) = 1 + 
^Var(Plf)+Vax(Pl^*) 
obtained with the a-parameter 

It is interesting that in principle any value of k will lead to creation of entangle- 
ment. The reason for this is our prior knowledge to the entangling pulse. Here 
the atoms are in the coherent state which is as well defined in terms of vari- 
ances as possible for separable states. We only need an "infinitesimal" extra 
knowledge about the spin state to go into the entangled regime. 
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It is interesting to see what happens to the conjugate variables Xai in the 
entangling process. This is governed by Eq. (25c). We do not perform mea- 
surements of the light operator Pj^- so all we know is that both X™- and 
are in the vacuum state. Hence Yav{X'^^) = (1 + k^)/2 and we preserve the 
minimum uncertainty relation Var(X^'J*)Var(P|^*) = 1/4. 

Entanglement Model With Decoherence 

Practically our spin states decohere between the light pulses and also in 
the presence of the light. We model this decoherence naively by putting the 
entire effect between the two pulses, i.e. we assume there is no decoherence 
in presence of the light but a larger decoherence between the pulses. We may 
then perform an analysis in complete analogy with the above with the only 
difference that P^ = PP^I + ^/l — P^Vp where is a vacuum operator 
admixed such that (5 = corresponds to a complete decay to the vacuum state 
and P = I corresponds to no decoherence. Completing the analysis we find 
the theoretical conditional variances 



1 J- n _ /92^ 2 

Var(A2|Ai) + Var(52|5i) = 1 + i^^—^, 

i + 

^ Var(Plf ) + Var(Pl^*) = ^\^y 
obtained with a-parameter 



(35) 



a = rr—:2- 



l + K^ 

In the limit j3 ^ 1 these results agree with (33) and (34). For /3 ^ we have 
a ^ (outcomes Ai and Bi are useless) and the variance approaches that of 
the vacuum state which is a separable state. 

4. Experimental Setup 

In this section we describe the details of the experimental setup, e.g. laser 
settings, pulse lengths, detection systems, etc. A picture of the experimental 
setup is shown in Fig. 1. In part (a) we see two cylindrical magnetic shields 
which each contain a paraffin coated vapour cell with cesium. The distance 
between the two cells is 35cm. In part (b) of the figure we show schematically 
the timing of laser pulses and the detection system setup. 

Laser Settings and Pulse Timing 

In one measurement cycle the first step is to create the coherent spin state 
of the atomic samples. To this end we have two diffraction grating stabilized 
diode lasers, one at the 894nm Dl transition 651/2 one at the 
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Figure 1. (a) A photographic view of the experimental setup. Atomic vapour cells are placed 
inside the cylindrical magnetic shields. The pumping beams are indicated with dashed arrows 
and the path of the entangling and verifying pulse is marked with the solid arrows, (b) A 
schematic view of the setup. The pulses reach a detection system measuring Sy{t). The pho- 
tocurrent is sent to a lock-in amplifier which singles out the sin(nt) and cos(f2t) parts. These 
are integrated and stored in a PC. The pulse sequence consists of (1) a 4ms pumping pulse, (2) 
a 2ms entangling pulse, a small delay r = 0.25ms, and (3) a 2ms verifying pulse. In addition 
to the pumping lasers and the entangling and verifying pulse, a laser beam is sent through each 
sample to measure the magnitude Jx of the macroscopic spins by Faraday rotation measure- 
ments. 



852nm D2 transition 6S1/2 6^3/2- We call these the optical pump and 
the repump, respectively. Both are sent through the first gas sample along the 
X-axis with a+ polarization, thus driving Am = +1 transitions only. In the 
second gas sample the polarization is (t_ . The main pumping is done by the 
optical pump laser, which drives the atoms towards the F = A, mp = 4 state 
(for (7+ -polarization). This state will be unaffected by the optical pumping 
laser (a dark state) because of the absence of an F = 5 state in the 6P1/2 
multiplet. In the pumping process some of the atoms will decay into the F = 3 
ground state, which is why we need the repumping laser from this state to the 
6P3/2 state to return them to the pumping cycle. Note that we can to some 
extent control the number of atoms in the F = 4 ground state with the power 
of the repump laser. With this optical pumping scheme we can obtain spin 
polarization above 99% (measured by methods similar to [19]). The pumping 
pulses last for 4ms and are shaped by Acousto Optical Modulators. 

The atomic samples are now prepared in the CSS with anti -parallel macro- 
scopic orientation. Next the off -resonant entangling pulse shaped by an Electro 
Optical Modulator of duration 2ms is sent through both cells to the Sy detec- 
tion. The difference signal is fed into a lock-in amplifier and the result is two 
numbers Ai and Bi corresponding to the integrated Jy^ + Jy2 and J'^^ + 

signals with an additional light contribution from the incoming S^^{t). The 
entangling (and verifying) pulses have power of P = 4.5mW and are blue de- 
tuned by 700MHz compared to the 6S'i/2, F = 4 ^ 6^3/2,^' = 5 transition. 
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They emerge from a Microlase Ti:sapphire laser which is pumped by an 8W 
Verdi laser. 

After the entangling pulse is a short r = 0.25ms delay before the verifying 
pulse is sent through the atoms. The entangling and verifying pulses have 
exactly the same shape and duration. The verifying pulse will again result in 
two numbers A2 and B2 stored in a PC. We now must predict the outcomes A2 
and B2 based on Ai and Bi as described in section 3. In order to calculate the 
statistics properly we repeat the measurement cycle 10.000 times. 

Measuring the Macroscopic Spin 

In addition to the entanglement creation and verification we also measure the 
macroscopic value Jx of the spin samples by sending linearly polarized light 
along the direction of optical pumping in both samples. This light experiences 
polarization rotation proportional to J^. As already noted we have for small 
angles 6 = Sy/2Sx- We also know Sy = aSxJx (for small angles, S^^ is zero 
classically) when light is propagating along the x-direction. Holding these 
together we find (this holds for all angles) 



^[rad] = — = -Y^;^- (37) 

Instead of the beam cross section A we need to use the effective cross section 
Agff = 6.0cm^ which is such that the volume of the vapour cell is F = ^eff • L 
where L is the length traversed by the laser beam (we call it "effective" since 
the vapour cell is not exactly box like). We only measure polarization rotation 
from atoms inside the beam cross section (which does not fill the whole vol- 
ume) and the equation must be scaled in order to count Jx for all atoms (we 

have J, = ji"beam.^^^/^)_ 

The DC Faraday rotation angle 6' is a practical handle on the macroscopic 
spin size Jx but we may in addition calculate a theoretical level for the pro- 
jection noise to shot noise ratio = PN/SN. Theoretically = a^JxSxT 
where again it is a question which cross section A to insert in a. The entangling 
and verifying pulses do also not fill the entire vapour cell volume. We assume 
that the measurement results are not depending much in the real cross section 
for the following reason. If A is small the interaction with the atoms inside the 
beam is stronger, but each atom will spend less time inside the beam leading 
to a correspondingly shorter effective interaction time T (atoms are at room 
temperature and move in and out of the beam). If this simple consideration has 
some validity we may assume that the light fills the whole cell volume and the 
correct cross section to insert is the effective area ^eff- Now this can be hold 
up against (37). Inserting A = 852nm, 7 = 5MHz and expressing the photon 
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DC Faraday rotation [deg] 

Figure 2. Atomic noise in units of shot noise as a function of atomic density (measured by DC 
Faraday rotation). Shot and electronics noise has been subtracted. Squares show the 1st pulse 
noise, circles the 2nd pulse noise. The linearity of these data is a finger print of the projection 
noise level which is then given by the solid line linear fit. Tip down triangles show the noise 
of the entangled states. Tip up triangles show the weight factor a. The two dashed curves 
trough triangles is the model described in Eqs. (35) and (36) with /3 = 0.65. The dash-dotted 
curve is the theoretically best for the triangles (/9 = 1). Note that the states created by the first 
pulse measurements are really entangled states (according to the criterion (31)) since the noise 
is clearly below the PN. We observe up to 25% noise reduction. Note, entanglement is observed 
for low densities also with ac^ < 1. 



tiux (j) = 2Sx 'iTi terms of probe power we may derive 

2 _ 18.6 • P[mW] • r[ms] • g[deg] 

'^theory- ^j^jj^j • (-^8) 

We remember the theory is a bit crude but we will see in Sec. 5 that the model 
holds pretty well. 

5. Experimental Results 

The experimental data is shown in Fig. 2, in the following we carefully 
explain the details of this graph. 

The Projection Noise Level 

The graph in Fig. 2 has on the abscissa the measured DC Faraday rotation 
angle which is proportional to J^- The squares are the variances Var(Ai) + 
Var(i?i) for the entangling pulses with shot noise (and electronics noise) sub- 
tracted. Also, the results are normalized to shot noise. The circles are the vari- 
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ances Var(yl2)+Var(Z?2) for the verifying pulses with the same normalization. 
According to Eq. (28) we have plotted the experimental ratio = PN/SN 
(the unity term in (28) is the subtracted shot noise). The linearity with Jx con- 
firms that we measure projection noise of atoms and not extra classical noise. 
We have roughly n^^^ = 0.10 ■ ^[dcg] which should be compared to the pre- 
diction (38) K^Y^^^^y = 0.24 • 0[deg]. The discrepancy is a little more than a 
factor of two but this is acceptable for our quite simple modeling. Note that 
the noise of the verifying pulse is the same as that of the entangling pulse. This 
is expected since we are performing a QND measurement. It is only when we 
remember the information given to us by the measurement results Ai and B\ 
that we can tell more about the state created by the entangling pulse. 

Conditional Variances and Entanglement 

The tip down triangles in Fig. 2 is the conditional variance Var(^2|^i) + 
\ai{B2\Bi) normalized to shot noise and with shot and electronics noise sub- 
tracted. According to (31) we thus plot k^{Ww{PXi) + Var(P|^*)). The fact 
that the points are lower than the straight Une (k^) is a direct indication that the 
entanglement criterion (27) is fulfilled. For the higher densities the reduction is 
25% but we note that entanglement is also observed for smaller densities with 

< 1. The latter was impossible with the older methods applied in [8]. The 
corresponding a-parameters from the minimization procedure (29) are plotted 
in Fig. 2 with tip up triangles. 

The expected entangled noise level in the ideal case is given by (33). This is 
drawn as the dash-dotted curve {n^ times 1/(1 + k^)). We see the conditional 
variance Ues higher than this curve and hence the entanglement is worse than 
expected. According to (34) we also would expect the a-parameters to lie on 
the same dash-dotted curve in the ideal case. It is clearly not the case, the 
experimental a-parameters are lower which indicates that the results Ai and 
Bi can not be trusted to as high a degree as expected. 

Let us try to apply the simple decoherence model given by Eqs. (35) and (36). 
Taking the decoherence parameter f3 = 0.65 we get the dashed lines in the fig- 
ure. These match nicely the experimental data. We conclude that the simple 
decoherence model has some truth in it and we must accept that the entangled 
state created can only be verified to be around "65% as good" as expected in 
an ideal world. 

Physical Decoherence Processes 

Above we quantify the observed decoherence with the /3-parameter. Here 
we comment on the physical grounds for the decoherence. A well known pa- 
rameter for describing the decay of the transverse spin components Jy and Jz 
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dt 



Ji, 
T2 



(39) 



where i = y,z. We have studied the T2-time extensively. A very good method 
for this is to create (by applying an RF-magnetic pulse) a displaced version 

of the coherent spin state with e.g. (JL) 7^ 0. This non-zero mean value can 



may be observed. 

Our experience tells us that power broadening by the laser pulses combined 
with light assisted atomic collisions play the important role in the decoherence 
processes. For high densities and high optical powers we may find T2 as low 
as 5ms. A fair guess for (3 is an exponential decay over a typical time scale of 
2ms (the time between the central parts of the entangling and verifying pulses). 
This yields /? ^ exp(— 2ms/5ms) « 0.67. This is not far from the observed 
P but we should say here that the 5ms is a typical value not directly measured 
in the case of the data given in Fig. 2. It is our experience that the observed 
decoherence in the entanglement experiments is stronger than that expected 
from the processes mentioned above. This is indeed true for shorter pulses. 
We believe that the atomic motion in and out of the laser beam combined with 
inhomogeneous hght/atom couphng due to the Doppler effect may also play a 
role, but we need further investigation to confirm this completely. 

Conclusion of Experimental Results 

To conclude the experimental section we emphasize that we have generated 
entangled states between distant atomic samples in the sense that each vapour 
cell sits in its own magnetic shield. The two shields can in principle be moved 
as far apart as is practical, our experiment was performed with a distance of 
35cm. In the future we hope to extend this distance further. The noise reduc- 
tion below the level set by separable states was measured up to 25%. This 
number is mainly Umited by power broadening and Ught assisted coIHsional 
relaxation of the atomic spins. The decoherence is successfully modeled by a 
single parameter (3. A much more elaborate theory on entanglement genera- 
tion in presence of decoherence and losses exists [17]. We should note, that 
the generated entangled states have random but known (based on yli and Bi) 
mean values. It is possible by applying RF-magnetic fields to shift the entan- 
gled states to having zero mean value while preserving the reduced variance. 
Experimental demonstration of such will be considered elsewhere. 




be detected by our standard /S^^-detection method and the decay following (39) 
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Figure 3. (a) Teleporting an unknown quantum state: first cells 1 and 2 are entangled. Then 
Alice sends a light pulse through cell 1 and the unknown quantum state in cell 3. The measure- 
ment results are communicated classically to Bob who by applying a displacement to his system 
based on the results of the two measurement recreates the unknown quantum state in his cell, 
(b) Entanglement swapping: the same procedure as in (a) but with two sets of entangled states 
initially. After a displacement based on the measurement results cells 2 and 4 are entangled 
even though they have never interacted directly. 



6. Perspectives 

We will now present a brief overview of some of the interesting applications 
in the field of quantum information of our reliable source of distant atomic 
entanglement. 

Atomic Teleportation. Quantum teleportation was first proposed in 1993 
[10] and the year after for the special case of continuous variables [20]. Tele- 
portation is extremely important since direct transport of physical states is of- 
ten hindered by exponential decoherence. With quantum teleportation the in- 
formation is cleanly separated into a classical part, which can be transmitted 
over arbitrary distances, and a quantum mechanical part, which only needs to 
interact locally. 

A proposal for spin state teleportation was given in [21]. Three atomic sam- 
ples are needed as shown in Fig. 3(a). Adjacent samples are oriented oppositely 
along the x-axis so that both collective measurements on cells 1 and 2 and on 
cells 1 and 3 will be regular entangling interactions as discussed earlier. Cells 
1 and 3 are located at Alice's site and cell 2 at Bob's site. The goal is now to 
teleport an unknown state in cell 3 onto cell 2. First cells 1 and 2 are entangled 
giving the measurement results {Ai, Bi). Next a pulse is sent through cells 1 
and 3 and the measurement results {A2, B2) are communicated classically to 
Bob. He can now by applying a rotation to his system based on the outcomes 
of the two measurements recreate the unknown state in his cell. This of course 
only works with perfect fidelity in the large interaction regime {k? » 1). 
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Teleporting an Entangled State: Entanglement Swapping. As described 
in [22] an entangled state can also be teleported using macroscopic samples of 
atoms. In this case Alice and Bob each have two samples (Fig. 3(b)). First 
each one of Ahce's samples are entangled with one of Bob's. Then a pulse of 
light is sent through Alice's two samples making it an entangled state. Alice 
now sends the result of a measurement on this last entangling pulse to Bob. 
Using this and the results of the two primary entangling pulses he can displace 
one of his states. This entangles his two samples without ever bringing the two 
into direct contact. Had Alice shifted one of her spin states prior to creating the 
entanglement between her two cells the exact same protocol would allow Bob 
to recover this shift in his cells. This allows for secret quantum communication. 

Light to Atom Teleportation: Quantum Memory. With an entangled 
source of atoms an unknown state of light can also be teleported onto this [11]. 
Assume we have two atomic samples with Jy\ + Jy2 = and Jzi + Jz2 = 0, 
i.e. a perfect EPR-entangled state. The protocol is simpler without rotating 
spins but also works with. As can be seen from Eq. (6) J, is mapped onto Sy 
when light propagating along the z axis is sent through an atomic sample. If 
light is sent through the first atomic sample and subsequently detected the mea- 
surement result can be fed back into Jz2 such that the original atomic variables 
exactly cancel. Sy has now been mapped perfectly onto Jz2- Another effect of 
the light pulse is seen from Eq. (8). Sz is mapped onto Jyi. If the transverse 
spins of sample 1 are rotated 90° and a new strong light pulse (k^ S> 1) is sent 
through this sample Jyf^ is measured. This result can be fed back onto Jy2 in 

such a way that the two original Jy's cancel. Sz is now stored in Jy2 and the 
teleportation is complete. This process could also be reversed so that an un- 
known atomic state is mapped onto a fight pulse via two EPR-entangled light 
beams. This shows that macroscopic samples of atoms offers a very feasible 
protocol for complete quantum memory. 

7. Summary and Conclusion 

We have presented the first experimental reafization of distant atomic en- 
tanglement in the sense that the two atomic systems are placed in separate 
environments, thus enabling entanglement between system separated by arbi- 
trary distances. Given the abundance of available quantum information proto- 
cols for this type of continuous variable atomic system the importance of this 
achievement is evident. Although quantum teleportation of atomic states has 
been achieved recently [23, 24] none of these approaches display directly scal- 
able distance between the unknown quantum state and the target system. Our 
system must therefore be considered an important candidate for achieving the 
long standing goal of high fidelity transfer of atomic states over great distances 
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upon which many of the proposed technical applications of entanglement and 
teleportation critically depend. The existence of several quantum information 
protocols for our physical system is based on the simplicity of the interaction 
between hght and atoms as expressed by Eqs. (25a-d). 
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